The Dark Matter Profile in the Galactic Center 
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We describe a quasi-equilibrium profile of dark matter particles in the inner parsec of the Galaxy, 
pdm oc r _3,/2 . This "mini-cusp" profile is caused by scattering with the dense stellar cluster around 
the supermassive black hole in Sgr A* and is independent of the initial conditions. The implications 
for detection of gamma rays from WIMP dark matter annihilation in the Galactic center are a mild 
enhancement of the flux, and a characteristic central feature in the angular distribution which could 
be detectable by high-resolution Atmospheric Cherenkov Telescopes. 

PACS numbers: 95.35.+d; 14.80.Ly; 98.35.Jk; 98.70.Rz 



The distribution of dark matter in the very center of 
our Galaxy is important for experimental searches for sig- 
natures of possible annihilation of supersymmetric dark 
matter particles. Such weakly interacting massive parti- 
cles (WIMPs) have remained favored candidate for the 
dark matter that represents most of the mass of the uni- 
verse ever since it was first proposed that the lightest 
supersymmetric partner particle, stable as a result of R- 
parity conservation, is the nonbaryonic dark matter 0. 
fn current standard supersymmetric theories these are 
expected to be neutralinos x El which naturally have 
the required cosmological density Q x ~ 0.25 0. 

Since dark matter makes a negligible contribution to 
the dynamical mass in the central parsec, its distribu- 
tion cannot be probed directly. Instead, it should be 
inferred considering all relevant physical processes in the 
Galactic center. Previous studies have reached different, 
often conflicting conclusions based on various initial as- 
sumptions and processes considered. All previous work 
assumed that dark matter particles are collisionlcss and 
therefore conserve their phase-space density. 

• On scales above ~ 1 kpc, most dissipationless sim- 
ulations of galaxy formation predict a power-law 
cusp in the dark matter density, pdm oc r -7 with 
7 = 1 - 1.5 0.0.0 or perhaps even 7 < 1 @. 

• In the vicinity of the Galactic center, the super- 
massive black hole with Mbh = 3.7 x 10 6 M Q dom- 
inates the mass in the inner r < 7'bh = 2 pc 
liol ITlj . If the central black hole grew adiabati- 
cally from a small seed, for example by accretion of 
gas, stars, and dark matter 0,0], the dark mat- 
ter cusp would be enhanced and would form a spike, 

oc r~ A with A = (9 - 2 7 )/(4 - 7) » 2.3 - 2.4 

• If instead the black hole appeared instantaneously, 
being brought in by mergers of progenitor halos, the 



enhancement is weaker and the spike has a slope 
A = 4/3 [l5j. A combination of this and the pre- 
vious effect, such that the mergers create a seed 
black hole that later grows by accretion, produces 
an intermediate slope. 

• Possible mergers of black holes in the centers of 
the progenitor halos 0, ^} may have the opposite 
effect on the dark matter distribution. Numerical 
simulations ^3 show that the kinetic heating of 
particles during the merger may reduce their den- 
sity to a very weak power-law, pdm oc r" 1 / 2 . How- 
ever, the simulations do not extend further than 
about 1 pc towards the center and thus cannot tell 
us about the very inner profile. 

• The same weak cusp, A = 1/2, results if the black 
hole grows away from the center of the dark matter 
distribution [III El]. 

Thus there is a considerable ambiguity in what the 
inner dark matter profile could be, all due to the uncer- 
tainty in the history of the central region of our Galaxy. 

Scattering off stars sets a quasi-equilibrium profile. — 
The above considerations assumed that the phase-space 
density of dark matter particles is conserved. However, in 
addition to the supermassive black hole, the Galactic cen- 
ter harbors a compact cluster of stars, with the density 
p* = 8 x 10 8 M Q pc" 3 HI in the inner 0.1" (0.004 pc at 
the distance of the Sun of 8 kpc) . These stars frequently 
scatter dark matter particles and cause the distribution 
function to evolve towards an equilibrium solution. Both 
stars and dark matter experience two-body relaxation. 

The idealized problem of a stellar distribution around 
a massive black hole in star clusters has been consid- 
ered in the past (cf. [23 f° r a review). Stars driven 
inward towards the black hole by two-body relaxation 
try to reach thermal equilibrium with the stars in the 
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core, but are unable to do so because of tidal disrup- 
tion or capture by the black hole. Unlike core collapse 
in self-gravitating star clusters, however, the density of 
inner stars does not grow toward infinity. A steady-state 
solution is possible where the energy released by removal 
of the most bound stars is transported outward by diffu- 
sion. Because there is no special scale in the problem, the 
quasi-equilibrium distribution function is a power-law of 
energy, f(E) oc E p , and the density is a power-law of 
radius, p oc r~ 3 ^ 2 ~ p j^,!^. The solution is unique and 
independent of the initial conditions. 

Within the sphere of influence of the central black hole 
rt,h, the distribution functions of both stars and dark 
matter are determined by two-body scattering and have 
the above power-law form. The scatterers in both cases 
are stars, but the distribution of dark matter differs in the 
exponent p from that of stars because of the vastly differ- 
ent masses of the two species. The evolution of the dark 
matter distribution f(E, t) in a two-component system of 
dark matter particles of mass m x and stars of mass m* 
can be described by a collisional Fokker-Planck equation 
(first derived in this form in [HI): 



dq df d 
~dE~dt " dE 
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where E = GMbh/r — \v 2 is the binding energy 
per unit mass, q{E) = (2 3 / 2 /3)7r 3 G 3 A/^ h ir 3/2 , A = 
167r 2 G 2 m 2 In A, and In A = In Mbh/m* ss 15 is the stan- 
dard Coulomb logarithm. The equilibrium distribution 
function of stars is f*(E*,t) oc E^ A , i.e. p = 1/4 [22| . 
For dark matter particles, however, the first term in the 
square brackets vanishes since the particle mass is negli- 
gible compared to stellar mass. An equilibrium solution 
with no energy flux requires df /dE = 0, or p = 0. The 
corresponding density profile is pdm °c r~ 3 / 2 . 

Dark matter particles cannot be tidally disrupted by 
the black hole, but they will be captured from the loss 
cylinder where their angular momentum per unit mass is 
less than J cap = AGM^/c (in which case their minimum 
distance to the black hole is within the last stable orbit). 
This would drive the outward flux of energy that needs to 
be balanced by the inward flux of particles, just as in the 
stellar case. As long as the relaxation time is shorter than 
the Hubble time, the quasi-equilibrium solution should be 
gradually achieved over a large range of energies. Note 
that this solution would be broken at the inner boundary 
where the distribution function would smoothly vanish 
due to the particle loss. The normalization of the profile 
may also slowly evolve (see for discussion). 

The observed stellar distribution around the black hole 
is roughly consistent with the equilibrium slope: p*(r) = 
1.2 x lO 6 (r/O.4pc)~ Q Af pc~ 3 with a = 2.0 ± 0.1 at 



r > 0.4 pc and a = 1.4 ± 0.1 at r < 0.4 pc 20]. The 
inner slope of the stellar profile is somewhat shallower 
than the predicted a = § +p = ?, but the overall profile 
is consistent with that of a relaxed cluster. The inner 
profile is well measured down to r = 0.004 pc from the 
center, where the local relaxation time is |25|: 



trel = 0.065 



G 2 m*p* In A 



2 x 10 9 yr, 



where v(r) « (GMbh/r) 1 ' 2 is the rms velocity of stars 
and dark matter, and m* « IMq. The density increase 
towards the black hole almost exactly balances the in- 
crease of the particle velocities, and the relaxation rate 
is roughly independent of radius. If the a = 1.4 stellar 
profile continues all the way to the gravitational radius of 
the black hole, r g = 2GM\>h/c 2 = 4x 10~ 7 pc, the relax- 
ation time is shorter than the Hubble time everywhere. 

The angular momentum of dark matter particles grows 
on the average as (GMbh?") 1 / 2 and the radius at which it 
equals J cap is £ C ap = 8r g w 3 x 10~ 6 pc. However, the 
extrapolated density at L cap is so high that dark mat- 
ter particles can annihilate faster then they are scattered 
in. For a typical value of the annihilation cross section 
(ff»f> ~ 10- 26 cm 3 s- 1 (c.f. Fig. 1), p = 100M Q pc" 3 
(see below), and m x ~ 100 GeV c -2 , the annihilation 
time equals the relaxation time at L w 10~ 5 pc. We 
therefore take this value of L as the smallest radius at 
which we can trust our dark matter profile. 

Note also that stellar-mass black holes resulting from 
stellar evolution of ordinary stars would sink towards the 
center from the inner 5 pc and would form a tight clus- 
ter of their own [2^|. This cluster scatters dark matter 
particles as well and additionally contributes to the re- 
laxation. 

Implications for dark matter searches. — The dark 
matter density in the central region of the Galaxy is thus 
given by 



Pdm{r) 



( I \-3/2 

Po(r/r hh ) 
Po (r/r hh y a 



L < r < r hh 



where L = 10~ 5 pc, and we expect that < a < 1.5. In 
order to calculate the flux of gamma rays due to WIMP 
annihilation near the Galactic center, we need to know 
the dark matter density po at radius rbh = 2 pc. This 
radius is at least two orders of magnitude smaller than 
the best currently available dissipationless simulations 
have reached. Unfortunately, the normalization of the 
dark matter profile in the central few pc is quite uncer- 
tain, since it is affected by various phenomena associated 
with the baryons which are the dominant mass compo- 
nent interior to the solar radius, d@. A starting assump- 
tion is that po is given by extrapolating inward an NFW 
profile 0j]. For example, the fit of corresponding to 
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p dm (d Q ) = 0.46 GeV 2 c" 4 cm" 3 with NFW scale radius 
parameter r s — 27 kpc, implies pdmfjbh) ~ 90 M Q pc~ 3 . 

Adiabatic c omp ression of the dark matter due to bary- 
omc infall is likely to result in an increased central 

density of dark matter. For example, model Al of the 
Galaxy in [i]) , including the effects of baryonic compres- 
sion, has pdm(100 pc) = 10 Mq pc~ 3 , with the density in- 
creasing inward as r~ 1A ; even if it scaled from 100 pc to 
rt,h only as this would give Pdm^bh)^ 560 M© pc -3 . 
On the other hand, in model Bl from 27], in which half 
of the angular momentum of the baryons is assumed to 
be transferred to the dark matter, pdm(100pc) is approx- 
imately 5 times lower than in model Al. In addition, if 
the black hole formed via mergers of the lower-mass black 
holes during the early stages of g alaxy formation, the 
density may be reduced further |18| . For definiteness, we 
take po — /0dm(7bh) = 100M Q pc -3 as our fiducial value, 
although it is likely that baryonic infall increased this. 

The flux at the earth from a region of angular radius 
ag = 0.05°, subtending solid angle AO = 7T(t|, centered 
on the black hole at the center of the Galaxy is then 



Flux = 



A 7 (a a nnV) I 

2d%m? x 



Here the number of gamma rays above the threshold en- 
ergy E t h of the detector is approximately [3(| 
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where x = E t h/m x c 2 and m x is the WIMP mass. As- 
suming E t h = 50 GeV and m x — 100 GeV c~ 2 , x = 0.5 
and N-y — 0.0087. Taking the distance of the sun from 
the center of the Galaxy to be c? Q = 8 kpc and as- 
suming a — 1 for definiteness, 



4 



"8"0 



2 

Pdn 



r 2 

Po d 2 



d e 



d.r 



n>h) 



= 0.73[24.4 + 5.2]pc 



V 100M Q pc- 



GeV' 



The minimum detectable WIMP annihilation cross sec- 
tion times velocity is then 



(cr a 



.1/2 
V bcg 



N 7 A cS tId- 2 



For an ACT with E th = 50 GeV, A bcg = (3.7 + 7.9) x 
10~ 6 cm -2 s _1 sr _1 A e g t AO, where the first term in the 
parenthesis is the electron-induced background and the 
second is the hadronic background. Since it is expected 
that the latter can be significantly reduced for the new 
generation of ACTs 0, |33| , we only include the electron- 
induced background. Then for a detection of statistical 
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FIG. 1: Minimum detectable annihilation cross section times 
velocity as a function of WIMP mass. The filled circles corre- 
spond to SUSY model WIMPs with Q x h 2 = 0.11 ± 0.01 
and the open circles correspond to SUSY models with Sl x h 
between la and 2a away from the central value. 



significance M s a (taking M s = 3 as our fiducial value), 



(Cann^mi!- 



3.2 x 10 -26 cm 3 s _1 
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The solid curve in Figure 1 shows this detection thresh- 
old as a function of the WIMP mass m x for our fiducial 
assumptions. The dependence on m x here comes entirely 
from m 2 /A 7 , but a more realistic calculation would take 
into account the increasing A c s for higher gamma ray 
energy. In addition, the significance of the detection 
would be increased by including the gamma rays from 
the full field of view of the ACT, which is 5° for the 
H.E.S.S. ACT array in Namibia 33], where the Galactic 
center passes nearly overhead. The dashed curve shows 
an order of magnitude improved sensitivity from baryonic 
compression increasing p a by a modest factor of 10 1 / 2 . 
Since at least that much baryonic compression is rather 
likely | 35| , it is quite plausible that annihilation from the 
Galactic center will be detected by new generation ACTs 
(H.E.S.S., CANGAROO III, MAGIC, and VERITAS) if 
the dark matter is actually WIMPs of mass m x > 100 
GeV c" 2 . 
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The best strategy for detecting WIMPs of lower mass 
probably involves using the all-sky map that will be pro- 
duced by the Gamma Ray Large Area Space Telescope 
(GLAST) satellite to look for the wide angle annihilation 
flux away from the Galactic center, which can be reliably 
calculated using dissipationless simulations || since it 
is hardly affected by uncertain physics such as baryonic 
compression. The recent reanalysis of the EGRET 
data toward the Galactic center using a gamma ray en- 
ergy dependent point spread function already constrains 
some SUSY models with low-mass WIMPs, subject to 
the uncertainties we have discussed about the dark mat- 
ter density toward the Galactic center. 

To summarize, if the WIMP mass m x > 100 Gev c -2 
and the central density is only a little higher than our 
fiducial value po — 100 M Q pc -3 due to baryonic com- 
pression, then the high angular resolution of the new 
ACTs such as the H.E.S.S. array should permit detec- 
tion of the pdm oc r~ 3 / 2 central cusp that is inevitably 
associated with the dense star cluster around the central 
black hole. 
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